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Abstract. The influence of nonperturbative fields on instantons in quantum chro- 
modynamics is studied. Effective action for instanton is derived in bilocal ap- 
proximation and it is demonstrated that stochastic background gluon fields are 
responsible for IR stabilization of instantons. It is shown that instanton size in 
QCD is of order of 0.25 fm. Comparison of obtained instanton size distribution 
with lattice data is made. 

1. Instantons were introduced in 1975 by Polyakov and coauthors [[[). These 
topologically nontrivial field configurations are essential for the solution of some 
problems of quantum chromodynamics. Instantons allow to explain anoma- 
lous breaking of U(1)a symmetry and the r/ mass spontaneous chiral 
symmetry breaking (SCSB) Taking into account instantons is of crucial 
importance for many phenomena of QCD (see |^| and references therein). 

At the same time, there is a number of serious problems in instanton physics. 
The first is IR inflation of instanton, i.e. divergence of integrals over instanton 
size p at big p. Second, quasiclassical instanton anti-instanton vacuum lacks 
confinement. 

The most popular model of instantons is the model of " instanton liquid" , 
which was phenomenologically formulated by Shuryak |6). It states that av- 
erage distance between pseudoparticles is R ~ 1 fm and their average size is 
p ~ 1/3 fm. Thus, p/R ~ 1/3 and vacuum consists of well separated, and 
therefore not very much deformed, instantons and anti-instantons. However, 
the mechanism for the suppression of large-size instantons in the ensemble of 
topologically non-trivial fields is still not understood. 

On the other hand, QCD vacuum contains not only quasiclassical instantons, 
but other nonperturbative fields as well. In this talk we will demonstrate that 
instanton can be stabilized in nonperturbative vacuum and exist as a stable 
topologically nontrivial field configuration against the background of stochastic 
nonperturbative fields, which are responsible for confinement, and will find 
quantitatively it's size. In this way, we will follow the analysis performed in 0. 

2. Standard euclidian action of gluodynamics has the form 

S[A] = ^J d i xti(F^ u [A]) = \J ^xF^[A]F^[A], (1) 

where F^-A] = <9 M v4„ — d v A^ — i[A^,Av] is the strength of gluon field and 
we use the Hermitian matrix form for gauge fields A^(x) — goA^(x)t a /2 and 
tr t a t b = S ab /2. We decompose A^ as A^ = A^ st + + a M , where A™ st is an 
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instanton-like field configuration with a unit topological charge Q T [j4 lnst ] = 1; 

is quantum field and is nonperturbative background field (with zero 
topological charge), which can be parametrized by gauge invariant nonlocal 
vacuum averages of gluon field strength]]. 

In general case effective action for instanton in NP vacuum takes the form 

Z = e -S e « [A'-] = f [£) ] /-SlA'^+B+a] \ 5 (2) 



where (...) implies averaging over background field B^. 

Integration over a M and B^ corresponds to averaging over fields that are 
responsible for the physics at different scales. Integration over a M takes into 
account perturbative gluons and describes phenomena at small distances, while 
averaging over B^ (formally interaction with gluon condensate) accounts for 
phenomena at scales of confinement radius. Therefore averaging factorizes 
Z = (Zi){Z 2 ) (see fj§ for details) and effective action appears to be sum of two 
terms, "perturbative" and "nonperturbative". Perturbative fluctuations were 
considered in |^|,|^] , and it was shown that in NP vacuum standard perturbation 
theory for instantons changes, which results in "freezing" of effective coupling 
constant. The perturbative part of effective instanton action in stochastic vac- 
uum S^A 1 " 15 '] was shown to be 

Here m* ~ 0.75m ++ ~ lGeV, where ++ is the lightest glueball. 
Thus, we have for effective instanton action 

S cS [A mst ] = Sl s [A mst ] + S^i [A inst ] (4) 
S e N ff p [A inst ] = _ i n (z 2 (B)) = - In (exp{-S[A inst + B] + S[A inst ]}) (5) 

3. We consider effect of NP fields on instanton, i.e. we evaluate (Z 2 ). In this 
work we make use of the method of vacuum correlators, introduced in works 
of Dosch and Simonov . NP vacuum of QCD is described in terms of gauge 
invariant vacuum averages of gluon fields (correlators) 

= (tr Gfj, lVl (x 1 )$(x 1 ,x 2 )G t _ l2l , 2 (x 2 )—G flnUn (x n )$(x n , x\)), 

( x \ . 

where G^ v is gluon field strength, and &{x, y) = Pexp i J B^dz^ is a paral- 

\ v J 

lei transporter, which ensures gauge invariance. In many cases bilocal approx- 
imation appears to be sufficient for qualitative and quantitative description 
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of various physical phenomena in QCD. Moreover, there are indications that 
corrections due to higher correlators are small jnj . 

Tensor structure of bilocal correlator follows from antisymmetry in Lorentz 
indices. It is parametrized by two functions D(x — y) and D(x — y): 



{g 2 GUx,x )GUy,x Q )) = (G 2 ) 



N? - 1 



< ^y^^pfW + ^^-{n p n p S ua + n v n a 5 w ) - (/a <-> f) j , (6) 

where G liV (x, x ) = $(x , x)G liV (x)^(x, x ), = z^/\z\ = (x - y)^/\x - y\ 
is the unit vector, (G 2 ) = {g 2 G p \ v G p \ u ) and, as it follows from normalization, 
D(0)+D(0) = 1. 

Bilocal correlator was measured on the lattice (see [jl2| and references therein) , 
and functions D{z) and D(z) were found to be exponentially decreasing D(z) — 
Aq exp(— z/Tg), D(z) — Aizexp(—z/Tg)/T g , where T g is the gluonic corre- 
lation length. Besides, according to lattice measurements A-y -C Aq [Ay ~ 
Aq/VS). Lattice data are presented in Table [TJ. SU{3) full stands for chromo- 
dynamics with 4 quarks, while SU(2) and SU(3) quenched mean pure SU{2) 
and SU{3) gluodynamics, respectively. 



Table 1: Lattice data [[12 


for bilocal 


correlator 


(G' 


; ), GeV 4 


T g , fm 


SU(2) quenched 


13 


0.16 


SU(3) quenched 


5.92 


0.22 


517(3) full 


0.87 


0.34 



To evaluate S^g we use the cluster expansion: 



(x 2 ) - (x)> 



(exp(x)) = exp I (a;) 



2! 



(7) 



In bilocal approximation we find S^f = >Sdia + h^dia, + 5 para + S\ + S%, where 



S, 



dia 



1 

1 



(fx ( tr 



{[A^B V ] - 



(8) 



Vra = -^4 / d 4 xd A y{tr(F^(x)G^(x))tr(F pa (y)G p<7 (y))) (9) 



51 = 7-1 d4xc ^y{ tT ( F ^ A ^ B ^x tT ( F pA A P, B . 

52 = 7 1 d4xd4y ^ tr (F ^ G ^ tr ( f p^ a pi B ^ 



(10) 

(11) 



We use notations Sdia (diamagnetic) and S'para (paramagnetic) for contribu- 
tions (H) and into interaction of instanton with background field. Next, 
S^f[ can be expressed through bilocal correlator (see Jt| for details), for in- 
stance 

1 1 

5dia = {< V1 N*- 1 / ^ J J W^( A ») 3 x ( 12 ) 

o o 

x [D((a - 0)x) + 2D((a - 0)x)] (13) 
4. We use standard form for instanton field configuration 

AT = 2t b R b ^J^±f f , (14) 

where matrix F&p ensures embedding of instanton into SU (N c ) group, b = 
1,2, . . .N% — 1; (3 = 1, 2, 3, r]" u are 't Hooft symbols. In singular gauge profile 
function f(z) satisfies boundary conditions /(0) = 1, /(oo) = and the classi- 
cal solution has the form f(z) = 1/(1 + z 2 ). Of course, real instanton profile 
in NP vacuum is different. The problem of asymptotic behavior of instanton 
solution far from the center \x\ 3> p c was studied in detail in Refs. [^3|-|l5|. 
Our numerical analysis shows that the value of p c is almost not affected by the 
asymptotic of classical instanton solution provided that condensate (G ) and 
correlation length T g have reasonable values. 

Numerical calculations show that S'dia(p) is a growing function of p, and it 
ensured IR stabilization of an instanton. Numerical results for instanton size 
distribution dn/d 4 zdp ~ exp(—S e g) and corresponding lattice data |l6) are 
presented in Fig. |l|. All graphs are normalized to the commonly accepted in- 
stanton density 1 fm~ 4 . Different lattice groups roughly agree on instantons 
size within a factor of two, e.g. p — 0.3 . . . 0.6 fm for SU(3) gluodynamics. 
There is no agreement at all concerning the density N/V. As a tendency, lat- 
tice studies give higher density and larger instantons than phenomenologically 
assumed. 

Using our model we find for (G 2 ) = 5.92 GeV 4 and T g = 0.22 fm that 
p c ~ 0.15 fm, which is less than phenomenological (~ 0.3 fm) and lattice 
results (in full QCD we find p c ~ 0.25 fm). However, we can present physical 
arguments to explain these deviations. Lattice calculations include cooling 
procedure, during which some lattice configurations of gluon field are discarded. 
This procedure can result in a change in gluon condensate (G 2 ), and thus 
instanton size distribution is calculated at a value of gluon condensate (G 2 ) coo i 
which is smaller than physical value (G 2 ). Therefore, lattice data for average 
instanton size p should be compared with our calculations for p c , performed 
at smaller values of (G 2 ). We show dependence of p c on (G 2 ) for several 



values of T g in Fig. ^. One can see that increase of (G 2 ) results in decrease of 
instanton size, and that effect is a result of nonlocal " diamagnetic" interaction 
of instanton with NP fields. 

We did not go beyond bilocal approximation in this work. As mentioned 
above, this approximation is good enough not only for qualitative, but also 
for quantitative description of some phenomena in nonperturbative QCD. In 
the problem under consideration there are two small parameters. These are 
l/g 2 (p c ) ~ 0.15... 0.25 and 1/N C , and there powers grow in each term of 
cluster expansion. Moreover, we made an estimate for the sum of leading 
terms in cluster expansion Jl7| , and found that IR stabilization stays intact 
(p c appears to be a little smaller). Thus, proposed model describes physics of 
single instanton stabilization in NP vacuum, not only qualitatively, but also 
quantitatively with rather good accuracy. 




Figure 1: Instanton density dn/cfizdp and lattice data Jlq |. SU(3) full (solid line), SU(2) 
quenched (dotted line) and SU (3) quenched (dashed line) 




Figure 2: Instanton size as a function of gluon condensate (N c = 3, Ne = 4) at T g = 0.2 fm 
(dotted line), T a = 0.3 fm (dashed line), T a = 0.34 fm (solid line) 
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